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PACS 61 . 43 . - j - Disordered solids 

PACS 62.40.+i - Anelasticity, internal friction, stress relaxation, and mechanical resonances 

PACS 63 . 20 . kp - Phonondefect interactions 

Abstract. - We calculate the scattering rates of phonons on two- level systems in disordered 
trigonal and hexagonal crystals. We apply a model in which the two-level system, characterized 
by a direction in space, is coupled to the strain field of the phonon via a tensor of coupling 
constants. The structure of the tensor of coupling constants is similar to the structure of the 
tensor of elastic stiffness constants, in the sense that they are determined by the same symmetry 
transformations. In this way, we emphasize the anisotropy of the interaction of elastic waves with 
the ensemble of two-level systems in disordered crystals. We also point to the fact that the ratio 
7;/7t has a much broader range of allowed values in disordered crystals than in isotropic solids. 



Introduction. — The "universality" of the glass-hke 
properties of amorphous sohds have been pointed out 
almost four decades ago [1]. Some of these properties 
are: heat conductivity which is almost independent of the 
chemical composition of the solid and proportional to T^ 
(where T is the temperature) , specific heat proportional to 
T, and a long-time heat release [1-3]. All these properties 
are described theoretically with reasonable quantitative 
accuracy by assuming that the amorphous solid contains 
dynamical defects that can be described at low temper- 
atures as an ensemble of two-level systems (TLS) [4,5]. 
Nevertheless, glass-like properties have been found also in 
disordered crystals [6-14] and quasicrystals, [15-17] only 
that in these materials they are not as universal as in 
amorphous solids and, even more, they exhibit anisotropy. 

The deep nature of the glass-like properties-and there- 
fore of the ensemble of the TLSs-remains elusive, despite 
of the long and intensive efforts invested into their study. 
This makes the study of disordered crystals especially in- 
teresting, since there, knowing the structure of the unit 
cell and its modifications due to disorder, we may know 
which are the tunneling entities and therefore we may 
have additional information about the TLSs. Moreover, 
the observed anisotropy of the glass-like properties, al- 
though unexplained, represents additional information for 
the theoretical description, which may help to improve the 
microscopic model. 



In general, the thermal properties of a dielectric glass 
are determined by the ensemble of TLSs, the phonon gas, 
and the interaction between them. In the standard tun- 
neling model (STM) the TLS is described in a basis that 
diagonalizcs the interaction Hamiltonian between the TLS 
and the phonon. In this basis, the Hamiltonian of the free 
TLS and the interaction Hamiltonian are 



H' 



TLS 



A -A 
-A -A 



and Hi = 



1 




-S 



(1) 



respectively. The interaction clement, S, is linear in the 
strain field of the phonon, [S] , namely 



6 = 2-fijSi, 



(2) 



where we assumed summation over the repeated sub- 
scripts. The symmetric second rank tensor [7] char- 
acterizes the TLS and its "deformability" under elas- 
tic strain. For the convenience of the calculations 
we work in the abbreviated subscript notations and 
we write [S] and [7] as the six-dimensional vec- 
tors S = (5ii, 522,5-33, 2^23, 25i3,25i2)* and 7 = 
(711, 722, 733, 723, 713, 712)*, respectively (where the super- 
script t denotes the transpose of a matrix or a vector). 
To go further and obtain deeper information about the 
TLSs, in rcf. [18] 7 was written as a product of two ten- 
sors: the first is a symmetric second rank tensor describ- 
ing the "free" TLS-we call it T in abbreviated subscript 
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notations and the second one is a fourth rank tensor, de- 
scribing the coupHng between T and S. In the abbreviated 
subscript notations, the fourth rank tensor is a 6 x 6 ma- 
trix which we shaU call [R] and represents the matrix of 
TLS-phonon coupling constants. So we write eq. (2) in 
matrix notations as [18] 



2T* • [R] ■ S. 



(3) 



The advantage of eq. (3) is that it separates the ten- 
sor T, which contains only the characteristics of the TLS, 
from the matrix of coupling constants, [i?], in which are 
embedded the characteristics of the interaction and has 
a general structure determined by the symmetries of the 
lattice [18, 19]. In general, the tensor T was taken in the 
simple form, T = (t^, ij, ^i, 2^2^37 2tit3, 2iii2)*, where ti, 
^2, and ^3 are the components of the unit vector t, which 
determines the "direction" of the TLS. This simplifica- 
tion will be used also in this paper. In the form (3), the 
symmetries of the lattice are imposed on the matrix [R] 
by coordinates transformations that leave the lattice in- 
variant, whereas the distribution over the elements of T 
is determined by the distribution over the "directions" , 
t. Through the properties of [R], this model predicts 
anisotropic glass properties of a crystal, even for an en- 
semble of TLSs isotropically oriented. 

In ref. [20] we applied this model to study the anisotropy 
of the glass properties in a disordered cubic crystal and we 
compared our calculations with the experimental results 
of Topp and coworkers [11, 13, 14]. Unfortunately the ex- 
perimental data published to date is not enough to check 
the model of ref. [18] or even to determine its parameters. 
From the available data, in ref. [18] we merely obtained a 
relation between these parameters, which should be con- 
firmed or not by future experiments. 

In this paper we extend our calculations to two other 
classes of crystal symmetries: trigonal 32 and hexago- 
nal. The former symmetry class corresponds to (neutron- 
irradiated) quartz and the latter to Na doped /3-AI2O3. 
Both materials show glass-like properties at low tempera- 
tures and strong anisotropy in the TLS-phonon coupling 
[7-9].^ 

It is known that in isotropic amorphous materials the 
coupling of TLSs with the phonon modes is described by 
the scalar coupling constants, 7; and 74, obtained by av- 
eraging the transition rates over the isotropic distribution 
of the TLS orientations. In this way, from very general 
considerations, one gets [18,21,22] 



(7//7t) > 4/3. 



(4) 



But in the model that we use here, this relation is affected 
by the symmetry of the lattice and therefore it does not 
necessary hold in a disordered crystal. This motivated us 
to discuss at the end of the next section the range of 7//7t 
for a crystal with cubic symmetry. 

Phonon scattering rates in trigonal and hexago- 
nal lattices. 



General considerations. The transition amplitude 
from a quantum state consisting of an unexcited TLS 
and nis^^ + 1 phonons of wavcvector k and polarization 
a {a — l,t), \n-i^cr + l,i), into the state of nto- phonons 
and excited TLS, |nk<TjT): is 

(nk.,t |i?i|nk. + l,i)=--V^^T*-[i?]-Ska (5) 



where e — y/A^ + A^ is the excitation energy of the TLS. 
Therefore the phonon scattering rate by a TLS in its 
ground state is 



rk.(t) 



2tt A^riko 



[i?]-Sk.|'<5(e-M- (6) 



The main characteristic of the TLS-phonon interaction is 
contained in the quantity Mk,(j(t) = T* • [R] ■ Skcr- As 
explained in the Introduction, the TLS-phonon interac- 
tion bear an intrinsic anisotropy through the matrix [i?], 
on which the symmetries of the lattice are imposed. To 
calculate the average scattering rate of a phonon by the 
ensemble of TLSs, we have to average over the distribu- 
tion of t. To reduce the number of degrees of freedom 
of the problem, in what follows we shall assume that t is 
isotropically oriented. 

Trigonal lattice. For a trigonal lattice of symmetry 
class 32 (the symmetry of quartz), the matrix [R] has the 
form [23] 



[R] = 



( rn 


ri2 


riz 


ri4 








ri2 


rii 


riz 


-ri4 








ri3 


ri3 


r33 











ru 


-rii 





r44 




















r44 


ri4 


\ 











7-14 


rii—ri2 
2 



\ 



/ 



(7) 



similar to that of the tensor of elastic stiffness constants, 
[c], with Cij replaced by nj [18, 19]. The system of co- 
ordinates that we use here is such that the z and x axes 
are the 3-fold and 2-fold rotational symmetry axes, re- 
spectively, while the y axis is perpendicular to both x 
and z. Solving the Christoffel equation we find that the 
crystal can sustain pure longitudinal waves propagating 
along the x and z axes, and pure transversal waves prop- 
agating along the y and z axes. The sound velocities of 
the longitudinal waves propagating in the x and z direc- 
tions are w^j — yjcwj p and v^^i = y/css/p, respectively, 
where p is the density of the material. The transversal 
waves propagating in the z direction have a sound velocity, 
Vz.t = \/C4a/p, independent of the polarisation direction. 
The pure transversal waves propagating in the x direction 
should be polarized only in the z direction and have a 
sound velocity va^t = \/{cii — ci2)/2/9-transversal waves 
polarized in other directions are not eigenvectors of the 
Christoffel equation. If we define the direction t by the two 
Euler angles 6 and (j), t = (sin0cos0, sinflsini/i, cos^)*, 
then for the longitudinal waves propagating in the x and 
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z directions we have 

Mki,i{0, (f) = ik[rii sin2(6') cos^{(t>) + ri2 sin2(6') sm^{(j)) 
+ri3 cos^ (6*) + ri4 sin(2^) sin{(f>)] (8a) 

and 

MkiA^, 0) = ifc[ri3 sin2(0) + r33 cos2(0)], (8b) 

respectively, whereas for the transversely polarized waves 
propagating in the y and z directions we have 

Mky,t{d,4>) = iks\n{9)[ruS\n{e)cos{24>) 

+2r44 cos(0) sin(0)] (8c) 

Mki:,t {0,4>) = 2ik sm{9) cos((/)) [r44 cos(6') 

+ri4sin(0)sin(</))), (8d) 

respectively. For the transversal wave propagating in the z 
direction we choose the polarization along the x axis-this 
choice becomes irrelevant after averaging over the direc- 
tions t, which we shall do next. 

The phonon absorption rates are calculated by averag- 
ing (6) over the ensemble of TLSs. If we denote by f{9, 0) 
the distribution over the angles of t{9, (/)), then we have 



we obtain the anisotropic values of the 7 parameters, 



and 



Potanhl2^^ 
T^„ = r^^ -n^a 



'kcr 



2h 



d9 

^0 



x|Mk.[t(0,0)]|2/(0,0) 
xnk.(|Mk,(t)|2). 



27tPq tanh 



2feBT 



(9) 



As mentioned before, we shall use the isotropic distribu- 
tion, f{9,(j)) — 1. Plugging eqs. (8) one by one into (9), 
we get the scattering rates 



'kx,l 






'ky,t 



' ki.t 



1 

15 



(Sr^i + 2riiri2 + 2ri3rii + 'ir^2 + 2ri2ri3 



-'irl^ 



4rf 



2TTPoN^nP 

n 



tanh 



2kBT 



8r^3 + 4ri3r33 + 3rl^ 27rPoJV^nfc^ 
15 ' h 



X tanh , , „ 
2kBT 

4(rf4 + r|4) 2TrPoN^nk^ 



4(4r 



15 

■2 4 
14 I 



2TrPoN^nP 



15 



tanh 
tanh 



2kBT 
e 



2kBT 



;(}Qa) 

(10b) 
(10c) 
^Od) 



where by N we denoted the normalization constant of the 
phonon {N — y/h/{2Vpujj) and by n the thermal popu- 
lation of the phonon mode (n = [exp(/3?ia;) — 1]~^); V is 
the volume of the solid. Comparing eqs. (10) with the 
standard formula for the phonon scattering rates. 



(STM)y 



2 2TTPoN^nk^ 
7,- tanh 



2kBT 



(11) 



Ikx.l 



lki.,1 



^ky.t 



2 

7fcz.t 



2{rl^ + rl2 + rl^) + (ni + n2 + r^f + Arf^ 



15 



iris + 4ri3r33 + 3rp 



33 



15 
4(r^4 + rji) 

15 
4(4r^4 + r|4) 

15 



(12a) 
(12b) 

(12c) 

(12d) 



Hexagonal lattice. The difference between the trigo- 
nal lattice of symmetry 32 and the hexagonal lattice is 
that ri4 and C14 are zero. This enhancement of symmetry 
allows propagation of pure longitudinal and transversal 
waves in all the three directions, x, y, and z. The sound 
velocities of the longitudinal waves in these three direc- 
tions are y^cu/p for x and y directions, and y^c^^/p for 
the z direction. For the transversal waves propagating in 
the X direction, the ones polarized in the y direction prop- 
agate with the velocity \/{cii — Ci2)/2p whereas the ones 
polarized in the z direction propagate with the velocity 
y^c^i/p. The transversal waves propagating in the y di- 
rection are similar to the ones propagating in the x direc- 
tion: the waves polarized in the x direction have a sound 
velocity of \/{cii — ci2)/2p, whereas the ones polarized in 
the z direction have a sound velocity of yc^^Jp. Finally, 
the transversal waves propagating in the z direction have 
all the same sound velocity, \fcAAjp- 

For the quantities M, we get 



Mfei,(^,0) 



Mfcx,y,t(^,< 
Muy,l{9,i 



Mky.i.ti^,' 
MkyM^,' 

Mki.i{9,^ 
MkiM9,, 



= ik[rii sin^(0) cos^((/)) + ri2 sin^(0) sin^((/)) 
+ri3cos2(0)], (13a) 



iA:sin2(6')sin(2</)) 



ni - ri2 



(13b) 
(13c) 



-^ri3cos2(0)], 
iA:sin2(6')sin(2</)) 



rii -ri2 



= «A:sin(20) cos(0)r44, 
= ik[rii siiA{9) siiA{(f>) + ri2 sin^(6') cos'^{(f>) 

(13d) 

(13e) 

(13f) 

(13g) 
(13h) 

(13i) 



iks[ii{29) sin((/i)r44, 
iA:(ri3sin2(6')+r33Cos2(6')), 
ik sin(26') cos(0)r44 , 
ik sin(20) sin((/))r44 , 



in obvious notations: the first subscript indicates the prop- 
agation direction while the second one is used only for 
transversal waves and denotes the direction of polariza- 
tion. We plug these formulae into eq. (9) to get 



'k.o- 



,B,^,2.PoNV ^^^^ 



(7- ) 



2kBT 



(14) 



where the superscript (H) stands for hexagonal and is used 
to make the difference between these quantities and the 
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ones calculated in the preceding subsection. As before, we 
get the 7 constants: 

( (H).2 _ 2(rfi + rfa + rj^) + (m + ri2 + r^^f 



15 



(7Si)^ 



/ (H) s2 _ 8r^3 + 4ri3r33 + 3r§3 

V Ik-v.ii.tl \ ikx.v.tJ 



15 



(15a) 
(15b) 
(15c) 

f^m )2 = (^(H) )2^(^(H)_ )2 = ( (H)_ ^2^l!i4 
V /fex,z,t^ \iky,z.t) V ffez,x,i/ \ ikz,y,tf -r r • 

(15d) 

We notice that the constants 7^^^;, 7^,^], 7^^|j and 7^. ^ 
are equal to jki,h lki.,h 7fcy,t and 7fcz,t, respectively (eqs. 
12) if in the last ones we set ru = 0. 

Range of ^i/^t ■ We notice also that in general the 
relation (4), valid for isotropic mediums, is not necessar- 
ily valid for crystals, which have lower symmetry. For 
the lattices studied above, the ratio 7^/7^ in any of the 
three directions has complicated expressions in terms of 
the components of [i?] . For the trigonal lattice \B\ has six 
independent components, whereas for the hexagonal lat- 
tice it has five. Therefore a discussion about the ranges of 
^f l^t f'3'' such symmetries would be too general to be of 
much use. 

The simplest lattice we can discuss is the cubic lattice; 
its [c] and [R\ matrices have only three independent com- 
ponents: Cii, C12, and C44 for [c] and rn, ri2, and r44 
for [i?]. The 7^ ^ and 7^ ^ (we use the superscript (c) to 
refer to the cubic lattice) constants have been calculated 
in Ref. [20] for longitudinal and transversal waves propa- 
gating in the (100), (110), and (111) crystallographic di- 
rections and for an isotropic distribution of TLS orien- 
tations. Using the results of Ref. [20], we calculate the 
ratios (7; /7t )^ for the waves propagating in the three 
directions mentioned above. Denoting Q = 7'i2/^ii and 
^ = r44/rii, we obtain 



\'t / (100) 
(c) J 

^^t / (110)4 

(c) J 
^^t / (110). 2 



3 + 4C + SC^ 



4^2 


2 + 6C + 7C^ + 4^2 


4^2 


2 -H 6C + 7(2 + 4^2 


(i-C)2 


5 + 20c + 20(2 + 16^2 



2[(1-C)2 + 2C2 



(16a) 
(16b) 
(16c) 
(16d) 



Note that in the (110) direction there are two transversal 
elastic waves, of reciprocally perpendicular polarization, 
propagating with different sound velocities. The isotropy 



condition for [i?] is ( + 2^ = 1, which sets the range of 
{li"^ /I'f^f to [4/3,00), as stated in eq. (4). If the lattice 
has lower symmetry, then ( -I- 2^ 7^ 1 and we introduce 
the parameter K to quantify the anisotropy, by imposing 
C, + 2K^ = 1; therefore K ~ 1 corresponds to the isotropic 
case. We calculate the dependence on K of the ranges of 
(7j /7j )2 in the three crystallographic directions of eqs. 
(16). Replacing ( by 1 — 2K^ into (16), we get 




15 

T 

15 

T ' 

15 
4K^ 



1 

e 



lOK 



-8K^ 
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lOK 


+ 7^2 


+ 1 
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10 


+ 7 + 


1 



(17a) 
(17b) 
(17c) 



45 



(111) 



4(2i^2 _^ 1) ^2 

4(5ii'2 + 1) 
^ 2K^ + 1 ■ 



30K 

210 + 



1 



(17d) 



Obviously, the condition K — 1 restores the isotropic 

equation for (ji/jt)'^ [18]. What is interesting to note is 

that all the equations (17) are quadratic in 1/^ and attain 

their minima at 

1 _ AK 

1 ■ ~ ~T' 
s mm '^ 

and we obtain the following constraints on (7;/7t)2 in the 
three propagation directions of the cubic crystal: 



(18) 
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^^t / (100) 
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(c)' 



^7t 
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7j 



(c)' 



(110)4 
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^7t 



(c) 
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(c) 



(110), 2 
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4K' 



3 ' 


i^2+3 


3 ' 


K^ + 3 


3i^2 ' 


4 



(111) 



2K^ + 1 ■ 



(19a) 
(19b) 
(19c) 
(19d) 



Now we can see that although for X = 1 all the conditions 
become identical, namely (7; /7t )2 > 4/3, for K =^ 1 
the lower limits of (7; /7t )2 vary differently. For ex- 
ample for K ':$> 1, the lower limits for (7/ /7t )?ioo) and 
(7/ /7t )nio) 1 become very big (ex K^), the lower limit 
of (7j /7t )nio) 2 converges to 1/3, whereas the lower 

limit of (7; /7t )/iii\ converges to zero. 

li K ^ 1, the situation is the other way around. The 
limit (19a) converges to zero, (19b) and (19d) to 1 and 
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4, respectively, whereas the limit value (19c) converges to [13 



infinity, like \/ K^ . Therefore in a cubic crystal, 7j can 
become smaller 7^ if the matrix [i?] deviates significantly 
from the isotropic condition. 

Conclusions. — We calculated the average phonon 
scattering rates on TLS in trigonal and hexagonal crys- 
tals, to emphasize the anisotropy imposed by the lattice 
symmetry. The parameters of the model may be obtained 
by measuring 7; and/or 74 in some crystallographic di- 
rections and this enables one to calculate the coupling of 
TLSs with phonons propagating in any other direction. 
The number of 7iS and 7tS needed, depends on the num- 
ber of independent parameters of the tensor of coupling 
constants, [i?], which is determined by the symmetry of 
the lattice. 

We showed that the allowed limits of the ratio 7^/7^ in 
crystals with different symmetries are different from the 
one imposed in isotropic materials, which is 7^/7^ > 4/3. 
In principle 7^/7^ in crystals may take any value. 

The calculations can be extended easily to disordered 
crystals of any symmetry. Moreover, although we used 
in our calculations an isotropic distribution over the TLS 
orientations, the comparison of our calculations with ex- 
perimental data would enable one to find if our assumption 
is true or not. If it is not true, one can determine, at least 
in principle, the distribution over the orientations of the 
TLSs. 
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